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Comment on "Low-dimensional Bose liquids: beyond the 
Gross-Pitaevskii approximation" 
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This is a comment on the work of Kolomeisky et al., Phys. Rev. Lett. 85, 1146 (2000). We 
point out that they are using the wrong form of the energy functional for one- dimensional fermions. 
We point out two possible forms of the energy functional, both of which can be derived from first 
principles but using different methods. One is obtained from the collective field theory method, 
while the other is derived from the extended Thomas-Fermi method. These two forms of the energy 
functional do not support the soliton solutions which are obtained by Kolomeisky et al. 
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The Gross-Pitaevskii (GP) mean-field theory replaces 
the bosonic field operator by a classical field <E>(r, t) |Q. 
This approach has been highly successful in describing a 
dilute gas of trapped bosonic atoms ||. However, some 
authors have pointed out that the usual potential en- 
ergy |<i>| 4 which arises from the zero-range pseudopoten- 
tial in three dimensions, needs to be modified in lower 
dimensions. Kolomeisky et al. || have proposed replac- 
ing the quartic term by a |<&| 6 term in one dimension. 
This can be motivated by considering a one-dimensional 
Bose gas in which particles interact pair- wise via a repul- 
sive 5-function potential. This is the Lieb-Liniger model 
which is exactly solvable [Q. The model has only one 
dimensionless parameter, namely, g = h 2 p/muo, where 
p — |$| 2 is the density, Uq is the strength of the S- function 
interaction, and m is the particle mass. In the limit of 
g ^ (called the dilute approximation), this system is 
equivalent to a gas of noninteracting fermions whose en- 
ergy density is given by 7r 2 7i 2 p 3 /6m. Kolomeisky et al. 
use this energy density to study a dilute Bose gas and to 



obtain a stationary soliton solution. 

Our main criticism of their work is that they are us- 
ing the wrong form of the second derivative terms in the 
energy functional for one-dimensional fermions, as indi- 
cated below. Since a dilute Bose gas with repulsive in- 
teractions in one dimension is equivalent to a system of 
noninteracting fermions, one way to proceed is to use the 
collective field theory (CFT) method to derive the Hamil- 
tonian in terms of the density variables. This was found 
long ago H . In the absence of an external potential, the 
energy density is given by 



H 



2m dx 3 



(1) 



where the density and phase fields p and are related to 
the order parameter field by $ = ^fpe %b '. The equations 
of motion for this system are easily obtained and 
these do not support any solutions in which p is time- 
independent and inhomogeneous. This seems to contra- 
dict the fact that a system of noninteracting fermions 
can have a time-independent and inhomogeneous den- 
sity profile 0, for instance, in the vicinity of a hard 



1 



wall. The exact density profile in both those cases has 
damped oscillations whose wavelength is of the order of 
the interparticle separation. This shows a limitation of 
the CFT method; the CFT Hamiltonian given in Eq. 
(|l|) gives the correct description of a system of noninter- 
acting fermions only for long-wavelength density fluctu- 
ations ||. A different derivation of the kinetic energy 
functional for nonintcracting fermions placed in an ex- 
ternal potential is based on the extended Thomas-Fermi 
method ||. In one dimension, to order h 2 , this gives the 
inhomogeneous terms in the kinetic energy density to be 
fi 2 /2m(-p' 2 jYlp + p" /3). This is different, in both sign 
and magnitude, from the term used by Kolomcisky et 
al, and has no solitonic solutions. In contrast to this, 
the Hamiltonian of Kolomeisky et al. does support sta- 
tionary solutions with inhomogeneous densities such as 
solitons. However, this does not by itself justify the addi- 
tion of the term h 2 p' 2 /&mp to the Hamiltonian (Q) which 
describes noninteracting fermions. 

Our second comment on the work of Kolomeisky et 
al. is that the dilute approximation is not obtained in 
any of the Bose condensates studied so far In the 
experimental systems, the equivalent of the parameter g 
is Na/aHOi where a is the scattering length, and ano 
is the harmonic confinement length. This equivalence 
follows because the particle density p in the center of 
the trap is of the order of N/clho, while the scattering 
length a is proportional to the strength of the pseudopo- 
tential in three dimensions and is therefore analogous to 
uq in the one-dimensional problem. Ref. j^j states that 
Na/dHO typically goes from a number of order 1 to sev- 
eral thousands, because a/auo is usually of the order of 
10~ 3 while N typically goes from 10 3 to 10 6 . Hence the 
experimental systems cannot be considered to be dilute 
in the Lieb-Liniger sense, and the mapping from interact- 
ing bosons to nonintcracting fermions in one dimension 
is not valid for such systems. 
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